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Qh! Abstract 

in 

We analytically calculate some thermodynamic quantities of an ideal <?-on 
gas obeying generalized exclusion statistics. We show that the specific heat 

> 

■ of a <7-on gas (g ^ 0) vanishes linearly in any dimension as T — > when the 
particle number is conserved and exhibits an interesting dual symmetry that 

■ relates the particle-statistics at g to the hole-statistics at 1/g at low temper- 
ed ' 

atures. We derive the complete solution for the cluster coefficients b[(g) as 

a ; 

S a function of Haldane's statistical interaction g in D dimensions. We also 

find that the cluster coefficients bi(g) and the virial coefficients ai(g) are ex- 
O ■ actly mirror symmetric (7=odd) or antisymmetric (7=even) about g = 1/2. In 



two dimensions, we completely determine the closed forms about the cluster 



5— i 

and the virial coefficients of the generalized exclusion statistics, which exactly 
agree with the virial coefficients of an anyon gas of linear energies. We show 
that the <?-on gas with zero chemical potential shows thermodynamic proper- 
ties similar to the photon statistics. We discuss some physical implications of 
our results. 
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I. INTRODUCTION 



Recently, there has been extensive interest in generalized exclusion statistics (GES) [pHUfl 
initiated by Haldane. It has been realized |[J,[|,[7|-[l2| that there are several models obeying 



GES in which interparticle interactions can be regarded as purely statistical interactions by 
generalized exclusion principle. GES defined by Haldane generalizes the Pauli's exclusion 
principle through the linear differential relation [|TJ 

M a = -Y,9«p&Np, (1-1) 

where Ad a is the change of the available one-particle Hilbert-space dimension when the 
number of added particles amounts to ANp and a and (3 indicate different particle species. 
This definition of statistics is independent of the space dimension and obviously interpolates 
between boson (g a/3 = 0) and fermion (g a/3 = 5 a p) statistics continuously. Note g a p can have 
arbitrary nonnegative values. 



An earlier form of fractional statistics-anyons is related to braiding properties of 
particle trajectories, which is a peculiar property of two dimensions unlike GES. GES su- 
perficially seems to have little to do with the braid-group notion of 2D statistics. However, 
several people have shown []I]|2],[7]-£| that an anyon gas in the lowest Landau level (LLL) 
satisfies the GES given by the statistical interaction g = a (braiding statistics parameter). 
Moreover, Murthy and Shankar have argued || that anyons obey the GES by relating the 
statistical interaction g to the high-temperature limit of the second virial coefficient (VC) 
of the anyon gas which shows a nontrivial dependence on 2D braid statistics [plj . 



As discussed by Nayak and Wilczek || , free anyons are not ideal g-ons but interacting g- 
ons. It will be an interesting question whether the ideal GES provides good approximations 



of the anyon model in which the Chern-Simons flux is attached to charged particles ||15|| . In 
this respect, it will be shown that ideal g-on statistics provides a good leading approximation 
of anyon statistics, i.e., the anyon statistics has a correct limit of GES, in the first-order 
perturbation-cluster expansion [fT7|,[n| and in the linear energies in a regulating harmonic 
potential [[0| . 



In this paper, we analytically calculate some thermodynamic quantities of an ideal g-on 
gas with g a p = g5 a p. For the single-particle energy given by e(p) = ap n , we verify that 
exclusion statistics satisfy the third law of thermodynamics as conjectured by Nayak and 
Wilczek ||, i.e., when the particle number is conserved, the specific heat of a g-on gas 
(g 7^ 0) vanishes linearly in any dimension as T — > and show that the specific heat at 
low temperature exhibits an interesting dual symmetry that relates the particle-statistics at 
g to the hole-statistics at 1/g. We derive the complete solution for the cluster coefficients 
bi(g) and thus determine the VCs ai(g) as a function of Haldane's statistical interaction g 
in D dimensions. We also demonstrate that the cluster coefficient (CC) bi(g) and the VC 
cii(g) are exactly mirror symmetric (/=odd) or antisymmetric (/=even) about g = 1/2, so 
that the semion (g = 1/2) which is the case of ID spinon as discussed by Haldane |1|J20| 
must exhibit a peculiar property, a 2 i{\) = &2i(|) = f° r ^ N. In the case of D = n, 
we completely determine the cluster and the virial coefficients of GES in closed forms, and 



the results exactly agree with the VCs of anyon gas obtained by Sen |T7| and Dasnieres de 
Veigy jB|. This result implies that a free anyon gas with a small a or with linear energies 
can be well described by an ideal g-on gas. We also show that the g-on gas with zero 
chemical potential shows thermodynamic properties similar to the photon statistics and all 
the thermodynamic quantities for a given volume depend only on the temperature and their 
dependences can be determined by simple dimensional arguments. Finally, we discuss some 
physical implications of our results. 



II. THERMODYNAMICS OF GENERALIZED EXCLUSION STATISTICS 

Consider a one-particle energy spectrum divided into a large number of energy cells e a , 
each of which contains G a independent levels and N a identical particles. In the limit of a 
discrete energy spectrum, G a would be the degeneracy factor. If d^ a is the dimension of 
the one-particle Hilbert space in the ath cell with the coordinates of N a — 1 particles held 
fixed, Haldane's definition of exclusion principle, Eq.(|l.l|), leads to 
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d Na = G a - g(N a - 1). 

Let W({N a }) be the number of configurations of the system corresponding to the set of 
occupation numbers {iV a }. Then an elementary combinatorial argument gives fHHH 

Under the constraints of fixed particle number and energy, 

A" = y n 

^ a (2.2) 
E = T e N 

J -' £—ia ^a 1 v a; 

the grand partition function Z is determined by Haldane-Wu state-counting rule fl2.ip as 
follows [^l[ 



Z=J2 W({N a })ex V {J2 N a(f^-^)/kT}, (2.3) 

{N a } a 



where \i and T are the Lagrange multipliers incorporating the constraints, Eq. ( |2.2| ), of fixed 
particle number and energy, respectively. 

The stationary condition of the grand partition function Z with respect to N a gives 
the statistical distribution n a = N a /G a of an ideal g-oia gas with chemical potential \i and 
temperature T as derived by Wu ||: 

n a = ~fA—, (2-4) 

W{Ca) +9 

where the function w(( a ) satisfies the functional equation 

w(C a ) 9 [l + wiU] 1 ' 9 = ( a = e^~^ kT . (2.5) 

As with usual statistics, in the thermodynamic limit, all the thermodynamic functions of 
GES can be evaluated around the set {A^} of the most probable occupation numbers. 

The Eqs.( ^.4p and ( |2.5p give correct solutions for the familiar Bose (g = 0) and Fermi 
(g = 1) distributions because w(( a ) = ( a — 1 for g = and w(( a ) = ( a for g = 1. Since 
the function Tg{w) = w 9 (l + w) l ~ 9 — ( a is a monotonically increasing function for w > 0, 



the algebraic equation ( |2.5|) has only one positive root for g ^ and T ^ and thus the 
statistics is unique for a given g. Also, n a decreases monotonically as g increases. In other 
wards, as g is increased from boson (g = 0) toward fermion (g = 1) statistics, the occupation 
number n a decreases for a fixed ( a . 

At T = 0, the ditribution ( |2.4| ) reveals a quite surprising phenomenon, such as formation 
of Fermi surface for g ^ : 

1/g, if e a < fi = e F , 
n a = < (2.6) 

^ 0, if e a > [i = e F , 

where the Fermi surface e a = €p is determined by the condition J2e a <e F G a n a = N. This 
property has important effects, to be demonstrated later, on the low-temperature thermo- 
dynamics of the g-on gas. Since the statistics at g = is the Bose- Einstein statistics and is 
well-known, we will not present a detailed analysis in this paper. 

As observed by Nayak and Wilczek and Rajagopal j|, Eq. (|2.5| ) displays a duality 
property that relates the particle-statistics at g to the hole-statistics at 1/g with a rescaled 

w((J;-) = l/w(( a ,g). (2.7) 
9 



In terms of the distribution n a in Eq.( [2.4|) , the thermodynamic potential = —kThiZ 
is given by 

n = -PV = -kT y G a In 1 + (1 ~^ n ° . (2.8) 
a I- gn a 

^From the Eq. ( |2.5|) , one can see the following properties on w(( a ): 

If e a -> oo, then w(C*) -> Ca = e (e ^ )/feT = -e £a/feT , (2.9) 

z 

and, 

e -e F / g kT if / o 

if e a -> and T -> 0, then w(C*) -> = <( (2.10) 



if = 0, 



where the fugacity z = e M//feT . 

In our calculations, we will use periodic boundary conditions over the system's volume 
V = L D ; thus, the D-dimensional momentum eigenvalue p of single-particle is given by 

2irh 
p=— n 

where the D-dimensional vector n takes only integer values. Then, in thermodynamic limit 
as V —>■ oo, a sum over p can be replaced by an integration: 

A state of an ideal <?-on gas can be specified by specifying a set of occupation numbers 
{n p } so defined that there are n p particles having the momentum p in the state under 
consideration. Of course, as pointed out by Nayak and Wilczek we must take the cell 
size of the single-particle energy spectrum to be sufficiently large, i.e. G a ^> 1, in order to 
safely ignore the specific complication of GES coming from negative probabilities. Then the 
occupation number N a of the ath cell will be the sum of n p over all the levels in the ath 
cell. With this in mind, the total energy E and the total number of particles N in Eq. fl2.2| ) 
can be represented as the sum over the momentum p: 

AT = 5> P (4), (2.11) 
p 

E = J2^M (2.12) 
p 

with n p {g) = N a /G a determined by Eqs. (|2.4| ) and ( |2.5| ). In the following, we will denote the 
statistical distribution n p (g) at statistics g as n g for simplicity. We assume the single-particle 
energy of an ideal g-on gas is given by 

e p = ap n , p=\p\. (2.13) 

Using spherical coordinates in D dimensions, the pressure P given by Eq. fl2.8| ) can be 
represented as the following integration in terms of the statistical distribution n g at g: 
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kT 



27r°/ 2 a-D/" roo , £-1 , l+(l-g)n g 
(27ra) D r(D/2)n JO Ufc fc ™ m i_ g „ 9 ; 

. kT ^ de £ D/ 



~ (2TTh) D r(D/2)D n '- L Jo u ' c c n g-> 
where the last step is obtained through a partial integration and Eq. 
the particle density p = N/V is given by 

r ~ (2nh) D r{D/2)n JO Ut t " n 9 

z dP 
kT dz ' 



(2.14) 
In the same way, 



(2.15) 



where the second step is obtained through a partial integration and by using the relation 



^yAlg _ ^drig 
de dz 



Equations ( 2.14 ) and Q2.15 ) can be written in the following standard form [21 



where 



P _ 1 Ag) f \ 
kT ~ \ dJ ^ [Z) 

A n 



| poo ^ 



(2.16) 
(2.17) 



and the thermal wavelength X D 



n 



(2wh) D q- r{D/2) 
2A{kT)% T tD/ny 



(2.18) 

By the Eq. (|2.17|) , one can determine 



the Fermi energy ep in Eq. (|2.6|) in terms of the density p: 



(2.19) 



^From the definition in the Eq.( |2.18D , we obtain the relation 



/From Eq.flTT^), the energy E can be also expressed as 



(2.20) 



E_DkT {g) 

which obviously exhibits the statistics-independent relation 



(2.21) 



nE 
DV 



(2.22) 



According to the fundamental principles of statistical mechanics, the entropy S of the 
system is given by 

where the second expression is obtained by using Eq. (|2.22|) . The density fluctuations in 
the grand partition function Z can be expressed as 

< N 2 > - < N > 2 = z—z—lnZ = —z—z—P 

oz oz kT oz oz 

= NkTpn T , (2.24) 



where the isothermal compressibility kt = —y(^)r = j(§p)t and Eq. fl2.15Q is used. It 



can be easily shown that the isothermal compressibility in the case of \x ^ is given by 

- = wjh'rAu- (2 - 25) 

Equation ( |2.24| ) shows that the density fluctuations of a g-on gas are vanishingly small in 
the thermodynamic limit, provided that is finite. 



III. GENERAL PROPERTIES OF GENERALIZED EXCLUSION STATISTICS 



A. General cases 



First, we will calculate the thermodynamic quantities, such as the total energy E, the 
specific heat Cy = (dE/dT)y at constant volume, the entropy S, and the isothermal com- 
pressibility Kt, for the single-particle energies given by e(p) = ap n at low temperature and 
high density, i.e., X D p 3> 1. 

At the low-T and high-p limit (g ^ 0), the Eq.( 2.18 ) can be evaluated by the so-called 
Sommerfeld expansion method |2T| after partial integrations with respect to y: 

(A 



fl 9 \z) 



II 



kT 

—yq, 
fi 



(3.1) 



where the numerical numbers C 9 at statistics g is defined by 



q 



°° dxx j^R^l 7 x = (e-fi)/kT 
dx 



(3.2) 



and Cq = — . Then Eqs.( |2.l7D and fl2.21|) are given by 



N = 1 n u D/n Y 

V \ D (kT)^T(D/n)D^ 



( D/n^ 



E In 1+M . , 

V " ~ \ D {kT)%T(D/n)D + n^ 



V 3 ) 
( 



kT 

{—yq, 



V J J 



kT 

{—yq. 



(3.3) 



(3.4) 



Let us rewrite Eqs.(|3.3|) and (|3.4|) in the following forms 



D/n _ D/r 

tp — fi 



3 = 1 V 



I D/n^ 



q 



e = E (^y+% 



kT 

j=l p 



3 


/ 






( 


1 + 








n 


q 


V 


3 


J 





(3.5) 



(3.6) 



where the Fermi energy ep is given by Eq.( |2.19|) and Eq = j^^NeF which is the ground 
state energy of the g-on gas at the given density. 

Here, we will convert the x-integration into a w-integration using Eqs.( [2.4|) and (|2.5|) : 



q 



dx x- 



dw 



+ / dw 



dw 



dx 




1 




(w + 


g) 2 


1 




(w + 


g) 2 


1 





{g\nw + (1 -g)\n(l + w)} j 
{g\nw + (1 - g) ln(l + w)} j 
{g\nw + {l- g) ln(l + w)} j + 0{e~ eF/9kT ). 



(3.7) 



(w + g) 2 

Since Wq ~ e ~ e F/gkT ^ ^ e seconc i integration can be certainly neglected, which is the 

reason we use the partial-integrated forms for N/V and E. The dual symmetry given by 

1 

Eq.( |2.7|) relates the coefficients C| at statistics g to the coefficients C- at ^ up to 0(e~ tF ^ 9kT ) 
as follows: 
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cf = (-i)v'- 2 c/ 



(3.8) 



which can be proved from Eg. ( |3.7|) by the substitution, w — > —. Note that only fermions 
(g = 1) have the exact particle-hole dual symmetry: Cj = for all odd j |21 |. After a partial 
integration, CJ (j > 1) becomes 

1 



= -^'-Hinwru =0 -j 



Using the integral formulas 



dw 



w(l + w) 



{g In w + (1 — g) ln(l + w)} 



j'-i 



In w 



u>(l + w) 2 
ln(l + ty) 



(ln r 

C(2) 



7T 



we can easily calculate Cf and C|. The results are 



dw 



ln(l + w) 
w(l + w) 



(3.9) 



(3.10) 
(3.11) 



Cf = 0, 
/~ig _ it 2 

The higher coefficients of can be obtained numerically 



(3.12) 



Ci 



C| 



-3.6062, Cf = 7.2123, C| 



-29.2227, Cf 



-4.8082, Cf = 14.4247, C| = -26.6973, Cf 



-5.4093, C 3 4 = 21.637, C 4 3 



-116.891 
-240.276 
-415.612, etc. 



(3.13) 



-25.9758, C% 

The above results clearly display the dual symmetry in Eq. (|3.8|) . 

The chemical potential /i, the total energy E, and the specific heat Cy are given by 

ix 2 ,D ^,kT 



/2 = e F 

E = Eq 



l-g 



6 n 

2 



1)(— ) 2 + 0(T 3 ) 



l+^(l + -)(-) 2 + 0(T 3 
6 n €f 



C v /k = g^(l + -)^kT + 0(T 2 ). 
3 n e F 



(3.14) 
(3.15) 
(3.16) 



Note that, at low temperatures, when the temperature increases, the chemical potential \i 
increases in the case of D/n < 1, remains stable in D/n = 1, and decreases in the case of 
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D/n > 1. This is consistent with the well-known results for a fermion gas when n = 2 [f2~2 . 



Using Eqs. ( |3.14| ) and ( |3.15| ), one can easily calculate the entropy 5* defined by Eq. Q2.23l ): 



S 1 ,D + n„ 

— = — E- uN) 

k kT D ^ ' 

= + ~)^kT + 0(T 2 ) (3.17) 

3 n 

~ CV/fc. 

When the particle number is conserved, the specific heat of a g-on gas (g ^ 0) vanishes 
linearly in any dimension as T — > 0. Therefore, as conjectured by Nayak and Wilczek j3|, 
GES obviously satisfies the third law of thermodynamics, which implies that the ground 
state of a many-body g-on system is nondegenerate. Our above results also agree with those 
of Nayak and Wilczek for the case of g = ~ and, in the case of g = 1 and D/n = 3/2, 



reproduce the well-known results for a fermion gas The linear dependences on T and g 
of the specific heat of a g-on gas can be roughly understood in the same way for an ordinary 
fermion gas. Of course, the key point is the existence of Fermi surface at low temperature. 
(The energy of particles excited over ep is of order gkT, and the number of excited particles 
is of order (kT/ep)N. Thus the total excitation energy amounts to g(kT/ep)NkT, so that 
Cy ~ g{kT / ep)Nk.) Note that the linear dependences on T and g of specific heat are 
satisfied even in D dimensions with e(p) = ap n . Also the isothermal compressibility k,t 



given by Eq. ( |2.25|) can be calculated by using Eqs. (|3.3j ) and (|3.14|) 



K T = g-T - + 1 47T/1 2 T 

n 2 



1-9 T - -1 — ) 2 + 0(T 3 
on e F 



(3.18) 



Note that the behavior of the leading term is linear in g and the correction of order T 2 
vanishes when D = n. 

Since the statistics at g = is a well-known bosonic one, we will present only the results 
for the statistics at g = for comparision: 

N 1 00 z l 1 z 

7 = F£iWrb (3 - 19) 

E DkT™ z l . 

y'-JoY.^- (3-20) 
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When T — ► or z — > 1, the above results exhibit the characteristic properties of Bose- 
Einstein condensation [[TL]], unlike the statistics at g ^ 0. The important point is that only 
the Bose-Einstein statistics (g = 0) exhibits Bose-Einstein condensation (D > 3). 

Now, we will investigate the thermodynamic propreties at high T and low p in D dimen- 
sions. This region corresponds to \ D p C 1 or z < 1. For this purpose, we will take the 
following expansion of the function fj) g \z) in Eq. (|2.18|) : 



fl 9 \z) = Y.bt\g)z l - 



(3.21) 



l=i 



According to the statistical relation ( 2.2(J| ), the coefficients tif\g) satisfy the following rela- 
tion 



bt l \g) = l-b?\g). 



(3.22) 



The following results can be obtained from Eq.( [3.2l|) : 



pi 1 oo 



i=i 



1 Mt ■ 1 



An A J=1 



(3.23) 
(3.24) 



(£+1) 

where &/(<?) = i> h n (g) is the Zth cluster coefficient. The virial expansion |21J of the equation 
of state is defined to be 



D „\l-l 



kT 



(3.25) 



1=1 



where a\{g) is called the Zth VC. By substituting Eq. ( |3.25| ) into Eqs. (|3.23| ) and ( |3.24| ) and 
requiring that the resulting equation to be satisfied for every z, we obtain 
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461 - 26 3 
a 4 = -206^ + I8b 2 b 3 - 36 4 

(3.26) 

a 5 = 2(56^ - 72b 2 2 b 3 + 96§ + 16& 2 & 4 - 26 5 ) 
a 6 = -672^ + 11206^3 - 315b 2 b 2 3 - 2806|6 4 
+6O6364 + 506 2 &5 - 56 6 



To get the Zth CCs bi(g), let us expand the distribution n{rf) as 

00 

n( V ) = J2 c i(9)v l , (3-27) 
where 77 = = ze~^. Then the coefficients b\ u \g) in Eq. (|3.21|) 



arc 



bf\g) = C M (3.28) 

It can be easily shown that for bosons q(0) = 1 and for fermions q(1) = (— 1)' +1 . Although 
the coefficients b\ f\g) or ci(g) for general g can be obtained by expanding the Eqs. fl2.4j ) and 
(2.5) in terms of we will show in the next section the CC b\ u \g) can be expressed by 



the closed form 



r<«) = I n (3.29) 



All the VCs, ((7) , can be obtained from the Eq. ([3.26|) because we have already known all 
the CCs, k(g), from Eq.(gjg). 



The first few terms of the pressure P and the particle density p in Eqs. ( |3.23|) and ( |3.24| ) 
respectively, are given by 



P = ~T,(Z 



kT 2g-l 2 9ff 2 -9^7 + 2 3 

TTT [Z n Z H „ Z H , 3.30 

D 2g - 1 2 % 2 - 9g + 2 3 

X p = z —z -\ o z H . (3.31) 

2 n 2 • 3 n 
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The virial expansion ( 3.25Q becomes 
PV = NkT 



4 D/n 3 l+£ 

In other words, the second VC and the third VC are 

2*7-1 



(3.32) 



02 
a 3 



(2g - l) 2 9g 2 -9g + 2 



3^f 



(3.33) 
(3.34) 



Note that 03, as well as 0,2, is statistics-dependent in D dimensions; 03 is statistics- 
independent only in 2D, or more accurately in the case of D/n = 1. Of course, when 
g = or 1, our above results reproduce the standard results of statistical mechanics. 

Using the results given by Eqs. ( |3.31| ) and ( |3.30| ), the entropy S in Eq.( |2.23| ) and the 
compressibility in Eq.( |2.25D in the Boltzmann limit are 



I ~ N [P±J1 - InA^p] 
k n 



+N[ 



D + n2g- 1 D 



A^p - ln(l + 



2(7-1 



In 4. r 9 ^- 9 ^ + 2 _ (2^7-1)% \ 2D p2] 
P [+i 2-3^ 4^ j (^) 2 J 



(3.35) 
(3.36) 



The first part in Eq.( |3.35| ) is just the Sackur- Tetrode equation pi| , and the second part 
is a correction due to quantum statistics. Note that this quantum correction has no linear 
dependence on p when D = n and that a sign change exists at g = 1/2. 

We will prove a very interesting mirror symmtery for the CC bi(g) and the VC ai(g) 
about g = 1/2. Let us start with the following relation derived by the Eqs. (|2.4|) and fl2.5|) : 



;i-<7n g ) s [l + (!-<?) 



1 1-9 



n a -e kT . 

z 



(3.37) 



Now, Eq.( j3.37|) can be rewritten as 



(3.38) 



where g = g — \- Oviously, one can see that the form of Eq.( |3.38D is invariant under the 
mirror transformation 
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(3.39) 



n 9 ~ n -9- 



The above mirror symmetry implies that, if n g (z) is a solution of Eq. (|3.38|) , then —n^ g (—z) is 
also a solution. Since the algebraic equation ( |3.38|) must have a unique nonnegative solution 
for a given g as remarked at the beginning of Sec. II, we can obtain the interesting symmetry 
property 

n_ g {-z) = -n- g {z). (3.40) 

Note that the above mirror symmetry about the distribution n g is generally satisfied for all 
temperatures and densities. 

The following mirror symmetry is the direct consequence of Eq. ( |3.40|) derived from 
Eq-CT: 



C2i+i{-g) = C2i+i(g) 

(3.41) 

C2i(-g) = -c 2 i(g)- 

The above mirror symmetry immediately leads to the mirror symmetry about the pressure 
P and the particle density p 

fD/ 9 n(- Z ) = -fD/n( Z ) 01 P(~ Z , H?) = ~ PM) 

' 1 (3.42) 

/fcjHO = -flU*) or P(-z,-g) = -P(z,g). 

n n 

^From the Eq.( |3.42j ), one can easily verify the following mirror symmtery about g = 1/2 of 
the CC bi(g) and the VC a^g) in Eqs. (pp-(pj): 



hi+i{-g) = hi+i{g), a 2 i+i{-g) = a 2 i +1 (g), 

(3.43) 

hi(-g) = -b 2 i(g), a 2 i{-g) = -a 2 i(g). 

In particular, semions with g = 1/2 exhibit a surprising property, a 2 i(\) = &2j(|) = for 
V/ G N. We can observe that our previous explicit results consistently show the mirror 
symmetries given by Eqs. (|3.40|) - (|3.43|) . 
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B. I)=ii 



In this section, we will consider the special case of D = n, which allows more analytic 
results. This includes the usual non-relativistic system with e p = p 2 /2m in two dimensions 
as a special case. In the case of D = n, we have an explicit solution |2j directly obtained 
from Eq.(|2T71) and Eq.(|2~5j): 



1 . l + wiz' 1 ) 

P = Tn m - 



ze 



-9>~ D P _|_ e -^°P = 1. 



(3.44) 



X D w{z~ 

This result will be very useful in our calculations of the CCs and the VCs. 

When D = n, there is no need to perform the Sommerfeld expansion used in the previous 
subsection. ^From Eq.( |2.21| ), the total energy of the system can be evaluated to be 

1 



v r°° 
E=— dw 



X D Jwiz- 1 ) w(l + Vj) 

V 

x° 



{fi + gkTlnw + (1 - g)kT\n(l + w)} 



, 1 + w(z x ) qkT l+w(z x " 
/iln , \ ' (In 



ln(l + w) 



) 2 + kT dw 

Jwiz- 1 ) W{1 + W) 



w{z L ) 2 ' w{z 

The above expression can be rewritten in a more elegant form by using Eq. fl3.44|) : 



(3.45) 



E 1 D 2 kT r* D P u 
y = ~ 2 9kTX p + ^ j o du—. 



(3.46) 



where e u — 1 = 1/w. 

First consider the low temperature and high density regions, i.e., X D p ^> 1. Since 



X D p 3> 1 and thus w(z 1 
evaluated up to 0{e~ x ° p 



<C 1, the integral in Eq. (|3.45|) or Eq.( |3.46| ) can be 



as 



D 



E = fiN- 



gkTX 



2 V 



-N 2 



n 2 kT V 
6 A^* 



+ 0(e 



(3.47) 



For g = 0, Eq.( p.47|) becomes 



^ 7T 2 kT V . \ D r 

E = —)fl +0{e 



(3.48) 



since p ~ —kTexp(—X D p) — > as T — > 0. For g ^ 0, Eq.( p.47| ) becomes 



E = -Ne F 
2 



3 €f 



(3.49) 
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Now the specific heat Cy at constant volume is given by 



S 1 NkT ' 0(e- xD ?), ]fg = 0, 



c v /k={ *wnD/*> P > ( 3. 50) 

g^^kT + 0(e- xD r), if^O. 
The above results, of course, coincide with those of the general case when D/n = 1. However, 
note that, as T — > 0, the chemical potential p remains constant, i.e., equals to ep and the 
specific heat has only a linear correction on T up to exponentially small corrections. This 



is a peculiar property of the energy spectrum with constant density of states [22|. 

We will now investigate a systematic way of calculating the CCs and VCs of GES in 
the Boltzmann limit using the explicit result in Eq. ( |3.44j ). In this region, Eq. ( |3.46| ) can 



be explicitly calculated in terms of the density p; thus, all the VCs can be determined. A 
straightforward calculation leads to 

^. p+ !£^ AV + ,|J^_ (AV -, (3.51) 

where B\ is the /th Bernoulli number (B 1 — — |, B 2 — |, 5 4 = — ^, etc.) and B 2 i+i = 
for I > 1 [^] . /,From the above equation, the VCs in Eq. ( p.25|) are found to be 



a x = i, 

2g-l 

02 = — ' (3.52) 

fl 2m - (2lTl)!> 

«2Z+2 = 0, for V/ > 1. 
Remarkably, the VCs, except the second VC, are (^-independent; moreover, all the even VCs, 
except (12(g), vanish. Of course, these properties of the VCs are consistent with the general 
properties depicted in the previous subsection. 

In order to calculate the CCs, bi(g), we will again use Eq.( ^.44j ). From Eq.( |3.44| ), we can 
deduce the following relation 

where p = X D p, 7 = 1 — g and -8/(7) are the Bernoulli polynomials ||23|| . Since z(p) can be 
expressed in closed form in terms of Eq.( [3.53| ), the series expansion p(z) in Eq. (|3.24j) is an 
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inversion of the power series z(p) p3f . Using the method of calculus of residues, the CCs 
can be expressed as 



l-b' 



(2) 



1 d 1 - 1 ffr 1 



l\dp 



2 = 



/! 



(3.54) 



where B«(i) are the generalized Bernoulli polynomials of order Z: 



t l e x 



(et-l) 1 

Thus, the CCs b, are 



n! d'- 
Z! 



^-(x — l)(x — 2) • • • (x — n), I >n. 



(3.55) 



»f>(9) = i- 2 n 



1 U 1 7/ - m 



m 



(3.56) 



Notice that the coefficients q(g) = bf\g) in Eq.( |3.28| ) are independent of order v. Hence, 
we can obtain a remarkably simple result about the coefficients b\ u \g) for general v from 
Eq.(p36D: 



1 ^ 7/ - m 



m=l 



(3.57) 



^From the above equation and Eq. (|3.26| ), the mirror symmetry of the CCs and the VCs in 
Eq. fl3.43p can also be directly proved. Interestingly, the similar results have been noticed in 
the Calogero-Sutherland model ||25|| , an anyon gas in the LLL [[F], and the anyon model with 



linear energies |19[ in a regulating harmonic potential. 

When D = n, we have completely determined all the CCs and the VCs and some lowest 
terms among them are 



bi = 






di = 


1 


b 2 = 


l-2. 9 
4 ' 




a 2 = 


2.9-1 
4 


h = 


32 ri m =i(i 


— —g), 


a 3 = 


1 

36 


w = 


4^ IIm=l(l 




a 4 = 





h = 


52 n m =i(i 




a 5 = 


1 

3600 



(3.58) 
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and a6 = 0, etc. The above results exactly coincide with the VCs (CCs), obtained by Sen 
fT7| and Dasnieres de Veigy |TjJ, for anyon gas. This result clearly shows that an anyon gas 
with a small a or with linear energies (g — 1 — a 2 ) can be well approximated by an ideal 
g-on gas. Thus the anyon gas confined in the LLL and the anyon gas with linear energies 
in a regulating harmonic potential can be regarded as the ideal g-on gas. Although the 
CCs show complicated statistics-dependences, the VCs are extremely simple and statistics- 
independent, except for the second VC. This implies that there are miraculous cancellations 
in the virial expansion. 



C. n = 

Let us consider the case in which the particle number is not conserved, p = 0. In this 
case, the number of particles N itself is not a given constant, but must be determined from 
the conditions of thermal equilibrium, (dQ/ '8N)t.v = p = The particle density p and 
the pressure P given by Eqs. ( |2.17| ) and (|2.16|) , respectively, become 

p = juf'o'il) = juI 9 d, 

(3.59) 

P _ i f(<?) \ _ l T g 

n 1 n 

where the definite integrals 1% and If , n are pure numbers. This shows that the pressure 

n n 

P and the particle density p depend only on the temperature T. Actually, for a given 
volume, all the thermodynamic quantities should be determined only by the temperature, 
and their temperature dependences should be obtainable by simple dimensional arguments. 
For example, the energy density E/V and the pressure are proportional to {kT) 1+ ^ and 
the particle density and the entropy density S/V to (&T)«. In the adiabatic expansion 
(or compression), the above results imply that the volume and temperature are related by 
VTn = constant and the pressure and volume by PV^ +1 = constant as in the photon statistics 
or an ordinary extreme relativistic gas. 

When D/n = 1, using the w-integration, the integrals if and if m Eq-( |3.59| ) can be 
expressed as 
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n 



9-1 



>1 



(3.60) 



where £ = ln(l + p) and p is a positive root of the algebraic equation, p 9 (l + p) 1 ^ 9 = 1. A 
numerical evaluation of some of the coefficients are as follows: 

l\ = 0.9624, If = 0.4812, If = 0.9870, If = 0.6580 

If = 1.1467, If = 0.3822, I| = 1.0785, If = 0.5664 (3.61) 
i? = 1.2891, if = 0.3223, l\ = 1.1400, I 2 4 = 0.5050, etc. 



IV. CONCLUSION 

We have analytically calculated some thermodynamic quantities of an ideal <?-on gas 
obeying generalized exclusion statistics. We have shown that the specific heat of a g-on gas 
(g 0) vanishes linearly in any dimension as T — > when the particle number is conserved 
and exhibits an interesting dual symmetry that relates the particle-statistics at g to the 
hole-statistics at l/g at low temperatures. At the low temperatures, the thermodynamic 
properties at g ^ are similar to those of fermion due to the existence of Fermi surface, 
unlike the case of g = (bosons) for which the Bose-Einstein condensation can occur. 

We have found a closed form for the CCs bi(g) as a function of Haldane's statistical 
interaction g in D dimensions. Unlike the low temperature case, high temperature behaviors 
are bifurcated at g — 1/2. The statistical interactions at g < 1/2 are attractive (bosonic) 
while those at g > 1/2 are repulsive (fermionic). The CCs bi(g) and the VCs ai(g) are 
exactly mirror symmetric (/=odd) or antisymmetric (Z=even) about g — 1/2. The case of 
g = 1/2 (semions) has a very peculiar property, d2i(\) = &2z(|) = ^ ^ or ^ e Spinon 
excitations in the Heisenberg spin chain with inverse- square exchange can be described by 
semionic statistics (g — 1/2) [fT|,p0|], and can also be related to half-filling Laughlin's boson 
fractional quantum Hall states J26|| . Thus the high temperature behaviors in these models 
may exhibit very interesting properties. 



20 



In the case of the energy spectrum with constant density of states, i.e. D = n, we have 
obtained full exact results and have calculated all the CCs and the VCs of a g-on gas and we 
have confirmed that they exactly agree with the virial coefficients of anyon gas, obtained by 
Dasnieres de Veigy |19[, for linear energies. Although this result implies that an anyon gas 



with linear energies can be well approximated by the ideal g-on gas and so the anyon gas of 
linear energy can be regarded as an ideal g-on gas, the free anyon gas cannot be generally 
regarded as an ideal g-on gas, e.g., in the case with a nonlinear energy spectrum or with 



strong coupling |L8| for which a statistical dependence of the third VC has been shown to 
appear in the second order correction. Even so, there may still remain the nonperturbative 
similarities, such as the mirror symmetry about semions, between anyon statistics and GES. 



Note that it was reported by Sen |27j that the third VC of anyons exhibits similar mirror 
symmetry about semions. It will be very interesting if the higher virial coefficients of an 
anyon gas turn out to be also mirror symmetric or antisymmetric about semions. 

In the case of \i = 0, the number of particles N is not a given constant as in a degenerate 
gas but the quantity being determined from the conditions of thermal equilibrium. As in the 
photon statistics, for a given volume, all the themodynamic quantities of a g-on gas depend 
only on the temperature and these dependences can be determined by simple dimensional 
arguments. 

In their paper ||, Nayak and Wilczek suggested that the electrons in Mott insulators 
behave as 2-ons and, thus, a Mott insulator has a Fermi surface of anomalous size and a 
specific heat of anomalous value. Our results show that the specific heat and the volume of 
the Fermi surface of 2-ons are twice as large as ordinary fermions in the units of the Fermi 



energy. Also, Ihm suggested |28| that the reduction of the number of allowed configurations 
of crystalline ice can be recognized as arising from the fractional exclusion due to the so-called 
ice rule and that this reduction can be described as the GES with the statistical interaction 
g = 0.867. However, it still remains open to construct an appropriate three-dimensionsional 
model describing the GES, unlike the one- dimensional soluble models nicely fitting with the 



notion of the GES such as Calogero-Sutherland model [lOR, antiferromagnetic spin chain 
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with inverse-square exchange [pO|j , an d the twisted Haldane-Shastry model fL2 |. Hence, it 
will be very interesting although very difficult to construct higher dimensional models nicely 
fitting with the notion of GES. 
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Note Added.: After the completion of the paper, we noticed Ref. [29], with which some parts 
of our main results overlap. 
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